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We investigate the signal from supernova relic neutrinos in future large scale observatories, such 
as MEMPHYS (UNO, Hyper-K), LENA and GLACIER, at present under study. We discuss that 
complementary information might be gained from the observation of supernova relic electron anti- 
neutrinos and neutrinos using the scattering on protons on one hand, and on nuclei such as oxygen, 
carbon or argon on the other hand. When determining the relic neutrino fluxes we also include, 
for the first time, the coupling of the neutrino magnetic moment to magnetic fields within the core- 
collapse supernova. We present numerical results on both the relic Ve and i^e fluxes and on the 
number of events for Ue+^^C, Ve+'^^Ar and Ue+p for various oscillation scenarios. The 

observation of supernova relic neutrinos might provide us with unique information on core-collapse 
supernova explosions, on the star formation history and on neutrino properties, that still remain 
unknown. 
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I. INTRODUCTION 



The observation of neutrinos from past core-collapse supernovae would represent an exciting and important discovery. 
Supernova relic neutrinos (SRN) bring crucial information on the explosion phenomenon, on the star formation rate 
and on neutrino properties. The understanding of how massive stars explode constitutes one of the main future 
challenges in astrophysics. While the various stages of the stellar evolution are quite well established, it still remains 
a mistery how, in these stars, the shock wave formed in the iron core, after collapse, manages to eject the mantle. 
Essentially all the gravitational energy is emitted in a short burst of neutrinos of all flavours. This luminosity curve 
follows closely the collapse, accretion phase and the cooling of the proto- neutron star left [l| . 

The first time such neutrinos have been observed is from the SN1987A located in the Large Magellanic Cloud, 
thanks to the 1MB Baksan Q and Kamiokande detectors J4|. The few electron anti- neutrino events measured 
have roughly confirmed the theoretical expectations, even though their energy and angular distributions still have 
unsolved problems. For the SRN fluxes, the best upper limits at 90 % C.L. of 6.8 x 10'^ i^e cm~^s~^ (25 MeV 
< < 50 MeV) and 1.2 De cm^^g-i > ^9.3 MeV) come from the LSD ^ and the Super- Kamiokande 

detectors The latter appears extremely close to most of theoretical predictions. Indeed it seems quite likely 
that the supernova relic neutrino background could be observed either by using an improved technology, e.g. the 
use of Gadolinium doped water Cerenkov detectors 0, or by scaling the present technology at megaton sizes, like in 
MEMPHYS 8], UNO ^9^, and Hyper-K [Tdj. These larg e-size observatories could provide us with the complementary 
information from a possible future (extra-)galactic explosion and from the SRN background. 

Many authors have estimated the expected SRN flux using different approaches to determine the supernova rate 



n , [13|, [IJ, [15|, [16|, [17[ (see [18| for a review). In fact, the latter can be deduced either from direct observations 
18 1 or from star formation rates (see e.g. (20j). The effects of neutrino masses and mixings on the SRN flux have 



been included recently (see e.g. [2l[), [24I). In [1^ a pragmatic approach has been used: to predict the SRN flux 
using available data, namely direct supernova rate measurements and supernova neutrino fluxes that give the best fit 
of the SN1987A events. Most of the available literature presents calculations for the De scattering on protons, which 
is the dominant detection channel in water Cerenkov (see e.g. [2^) or scintillator detectors [l3|- In particular, a 
study of what could be learned on the neutrino spectra and/or the core-collapse supernova rate at next generation 
water Cerenkov detectors is realized in [2^ . In Ref. [11] the potential of liquid argon detectors to observe supernova 
relic electron neutrinos is analyzed, while the specific case of i^e detection in SNO is considered in [13, [1^ • Ref. 28 1 
also obtains an upper limit on the relic I'e flux from the De upper limit in Super-Kamiokande. Finally, in 2^, 3C| 
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the authors have considered the possibility of using the SRN to gather information on more exotic unknown neutrino 
properties such as neutrino decay. 

The study presented in this paper has several new features compared to the available literature. First we study 
the possible signal from both the supernova relic Vg, and v,, in future large-scale observatories. We consider the 
three technologies under investigation, water Cerenkov, scintillators and liquid argon, having in mind the proposed 
MEMPHYS (UNO, Hyper-K), LENA and GLACIER detectors We also include for the first time in the neutrino 
propagation within the supernova the interaction with matter (the Mikheev-Smirnov-Wolfenstein effect as 
well as the coupling of the neutrino magnetic moment to the magnetic fields (the spin- flavour conversion phenomenon). 
Clearly our present understanding of supernovae explosion is still limited from several respects and, in particular, for 
the knowledge of the stellar magnetic fields and their evolution during the explosion. We think however that, even 
at the present stage, it is interesting to perform a first calculation including this process, to start investigating under 
which conditions some effects might be seen. We make here the assumption of a static supernova density profile, 
neglecting in this first work the changes induced by the shock wave. We also do not consider earth matter effects 
that have been discussed in To determine the SRN flux, we integrate the neutrino spectra from one supernova 
explosion over the redshift by convoluting it with the supernova explosion rate. We predict the expected event 
rates and analyse the effects of the presence of a neutrino magnetic moment both on the relic number fluxes and 
on the number of events. Note that we follow previous investigations for the choice of the neutrino fluxes at the 
neutrinosphere, of the supernova rate as well as for the discussion of the possible backgrounds in the detectors. The 
paper is organized as follows. The formalism is described in Section II. The results and conclusions are presented in 
Sections III and IV. 



II. FORMALISM 

A. Neutrino evolution in a core-collapse supernova with resonant spin-flavour conversion 

Neutrino evolution including coupling with matter and with stellar magnetic fields (the resonant spin-flavour conver- 
sion or RSF) was first considered as a possible solution of the solar neutrino deficit problem [s^jH^,!!^. We now know 
that this deficit is due to the presence of neutrino masses and mixings [37| . Even though the RSF might still give a 
sub-dominant contribution, it turns out to be extremely small [3S] . The situation could be quite different in the case of 
massive stars such as core-collapse supernovae (SNII), since strong magnetic fields are expected, particul arly close to 
the proto-neutron star surface. Many authors have investigated the RSF effects in SNII (see e.g. ^391 [iol. [ill [42l. |43| ) . 

The knowledge of the neutrino magnetic moment offers a serious challenge^. While the minimal extension of the 
Standard Model predicts a very small value 1], the observation of a large value would clearly indicate new physics. 
Several bounds have been obtained so far. Direct upper limits come from neutrino-electron scattering using neutrinos 
from reactors 37, 44, 45, 46], the best value being 9. x 10~^^ at 90% C.L. Based on the electroweak radiative 
corrections to the neutrino mass, a model- independent theoretical upper bound is derived in [47j . Important indirect 
limits are obtained from astrophysical observations such as red-giant cooling [4^ the SN1987A [s^], or Big-Bang 
nucleosynthesis (49j . One might hope to attain improved direct limits with running experiments at reactors, a very 
intense tritium source, low energy beta-beams [s^j or future astrophysical observations. 

The neutrino evolution equations for three flavours including the interaction with matter and the coupling to the 
magnetic fields is given by (?i = c = 1): 

where v [v) stand for the electron, muon and tau (anti-)ncutrinos. The hamiltonian includes the neutrino mixing 
and interaction with matter: 

, /O \ /Kc + Kc \ 

^0 = C/t + K'^ + 14c , (2) 

\Q ^mlj \ Kc + Kc/ 

where the first term corresponds to neutrino propagation in vacuum while the second includes V^^, the potential due to 
the charged-current neutrino scattering with electrons and y„c, the neutral-current neutrino scattering with nucleons. 



Note that since neutrinos are massive particles they certainly have a non-zero magnetic moment. 
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In Eq.([2]) U stands for the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix relating the neutrino flavor and mass 
basis. For anti-neutrinos, Hq is similar to Hq but with a minus sign in front of the second term. The off diagonal 
terms come from the neutrino interactions with the SNII magnetic field, B{r) and M is the magnetic moment matrix. 
Here we consider the case of Majorana neutrinos for which M is 



with the transition magnetic moments. (In [43| we have shown that neutrino flux modifications of less than 1 % 
are induced by the contribution of the neutrino magnetic moment to the neutrino scattering with matter. Here we 
do not consider these contributions.) From the RSF point of view Eq.ll]) the relevant quantity is the product BM. 
The observation of RSF effects from a future SNII explosion signal, or from supernova relic neutrinos, might provide 
us with further information on the value of ^i, and/or the strength and profile of the magnetic fields in SNII. 

Coupling with matter or with magnetic fields can induce important modifications to neutrino oscillations. While 
the former produce versus and v^., the latter can produce versus and resonant conversions. The 
MSW resonances are situated at two densities governed by (Amjgj^ia) and (Amf2,^i2) respectively. The first 
resonance is located at density of about 10^ g/crc? (MSW-/i), while the second one at about 1 g/cm^ (MSW-Z). 
The neutrino conversion at each resonance depends on its adiabaticity. While the MSW-Z resonance is adiabatic 
for typical supernova density profiles, the adiabaticity of MSW-/i depends on the still unknown mixing angle ^13 
[51I l52 |. Similarly, a resonant Vi versus Vj conversion (and vice- versa) can take place due to the presence of the 
neutrino magnetic moment. Three supplementary resonances are present in this case, their adiabaticy depending on 
the strength of the magnetic fields (see e.g. [l^l for a detailed discussion). 



Core-collapse supernovae are one of the most explosive phenomenon in our Universe. When the iron core of a massive 
star reaches its Chandrasekhar limit, an increased gravitational pull due to a contraction is no longer counterbalanced 
by the increase in the thermal pressure. The iron core collapses until it reaches nuclear densities and becomes very 
stiff. A shock wave is formed that moves outward ejecting the stellar mantle. How such ejection takes place remains 
still an unknown issue. Convection or extra energy deposit by neutrinos might offer the clue to this longstanding 
problem. In fact, most of the energy produced by the explosion is taken away by a short burst of neutrinos of all 
flavours. Their luminosity curve follows closely the different phases of the explosion. 

Two ingredients are crucial to predict the SRN differential flux, dF/dE,^, i.e. the neutrino flux from one supernova 
and the supernova rate : 



where TZsN{z,^l\,^m) is the supernova rate per comoving volume at redshift z and flm are cosmological pa- 
rameters, fl\ + Qm = 1 in the standard ACDM model), z^ax is the redshift at which star formation started and 
dn^/dE^{{l + z)E^) is the energy spectrum of emitted neutrinos leaving the supernova ^21]. 

The neutrino spectra dn^/dEi, at the neutrinosphere are quite well approximated by Fermi-Dirac [ij or power-law 
distributions [53^. Electron (anti)ncutrinos interact with matter via charged as well as neutral current, while muon and 
tau neutrinos via neutral current only. The matter being neutron-rich and the cross sections for anti-neutrinos weaker 
than for neutrinos, we expect the following neutrino hierarchy: {Ei,^) < (-Ep^) < (-Ejy^^p^) with x = ^,t. Supernova 
simulations predict typical temperatures of 10-13 MeV, 13-16 MeV and 16-23 MeV depending on the microphysics 
[5^. To determine the dn^/dEi, at the SNII surface, we solve the neutrino evolution Eqs.([l]|3]), where one of the 
im por tant unknowns is the stellar magnetic field. Here we consider that its shape follows a power-law, as done e.g. 
in liOl 




(3) 



B. Core-collapse supernova neutrinos and the supernova rate 




(4) 




(5) 



where n = 2 or 3 and tb = 1.4 x lO^cm, the iron-core surface. 
By including dt/dz and using 



nsN{z,iixM = (1 + z)-^'^^{i + nmz){i + z^) - nx{2z + Z^)nsN{z) 



(6) 
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FIG. 1: Cross sections as a function of neutrino energy for + p ,54] and Ue+^^C Ve+^^O [Sq ] and Ue+^'^Ai [57 



one gets that dF/dE^ is independent from the cosmological parameters: 



dF/dE^ 



1 

Hq Jo 



„ , , dn^((l + z)E^) ,^ , o/o , 
dEu 



(7) 



where we take Zmax equal to 4. For the supernova rate we consider the functional form from indirect observations 
accounting for the star formation rate 41 1 



UsNiz) = 1.22 X 10^2 X 0.3/^/16 



exp (3.4z) 
'exp(3.8z) + 45 



yr 



'^Mpc' 



(8) 



with /ig5 = i/o/65km s^^ Mpc The factor 1.22 x 10^^ is the fraction of SNII, deduced by assuming the Salpeter 
IMF (with a lower cutoff at 0.5 Mq and taking all stars with M > 8Mq as SNII). The normalization factor fi, 
takes into account the uncertainty on the local star formation rate density, which at z = is (0.5-2.9) x lO^^/iyoAf© 
year~^Mpc~^. The value = 1, that we will take in the following, is consistent with mildly dust-corrected UV data 
at low redshift [l^ and predicts a local star formation rate of 0.7x lO~^/i7oM0 year~^Mpc~'^. Note that according to 
[20| . the local value for the star formation rate is 1.6x lO~^/i7oM0 year~^Mpc~^ which means = 2.5. Here we do 
not perform a study of the sensitivity of the SRN flux to different choices of the TZsn{z) since this was performed in 
[2ll |. The different functionals available in the literature essentially have a very similar increasing behaviour for low 
redshifts (z < 1), well constrained by observations; while they might differ at high redshifts (z > 1). Note however 
that the latter contribute little to the signal in an observatory, if the detection threshold is larger than about 10 MeV. 

Once the relic fluxes are obtained, the event rates are calculated by multiplying the SRN flux by the corresponding 
cross sections (Figure [1]): 



NpdF/dE^a{E^)dE^ 



(9) 



Eq 



with Eq the threshold energy that depends on the detector technology, Np the number of target nucleons or nuclei, a 
the relevant cross section. 



III. RESULTS 



We have calculated the neutrino evolution in a core-collapse supernova Eqs.lUJS]) and then the relic supernova neutrino 
flux Eqs. (|4ll8)) . The osciUation parameters are fixed at Amfg = 8 x 10"^eV^, sin226'i2 = 0.86 and Am^g = 3 x 10"^eV^, 
sin^26'23 — 0.96 for the solar and atmospheric differences of the mass squares and mixings, respectively [33|. For the 
third still unknown neutrino mixing angle ^13, we take either the present upper limit sin^20i3 = 0.19 at 90 % C.L. 
(L) or a small value of sin'^26'i3 = 10"* (S) that might be attained at the future (third generation) long-baseline 
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FIG. 2: Density profile of tfie pre-supernova progenitor used in tfiis work |59i] . 



experiments [58| . Since the sign of the atmospheric mixing is unknown, we consider both the normal (NH) and 
inverted (IH) hierarchy, i.e. Am^g > or < 0. 

For the neutrino spectra at the neutrinosphere, we take Fermi-Dirac distributions as an example (Table [T|. We 
do not perform a detailed study of the sensitivity of the SRN flux predictions to different choices of the neutrino 
temperatures, from available supernova simulations, since this was done in [24j . The total luminosity is fixed at 
Lq = 3 X 10^^ ergs and the energy is assumed to be equally shared among all flavours (this is not necessarily true at 
early times). For the SNII matter density profile, we use a static 13.2 Mq progenitor [53| (Figure [2]). With the goal 
of also investigating the RSF effects on the SRN signal, we consider two different power law profiles Eq. ([5]) and the 
range of 0.001-1 values for the /iy-Bo product Eqs.((l][3]). These values correspond to taking at r = for example a 
magnetic field of Bq ~ 10^ {Bq — 10^^) Gauss and a transition magnetic moment of 10^^^ (10^^^) /is or Bq — 10^^ 
Gauss and 10~^^ /is respectively. In the following we will denote the largest value of 0.1 (1) for n = 2 {n — 3) by 
{fj-B)iarge- Notc that the radial profiles and the fJ-ijBo product are such that, while fj,ijB(rB) = 0.1, 1 for n = 2, 3; for 
both profiles fj,ijB{10rB) = 0.001 at the location of the RSF-/i resonance. However, the behaviour is different at the 
RSF-^ location. This has an impact on the expected rates, as we will see. 



Electron fraction Ye 








(S.J 


{E.J 




0.48 


3.5 MeV 


5.1 MeV 


7.0 MeV 


11 MeV 


16 MeV 


22 MeV 



TABLE I: Neutrino temperatures and average energies at the neutrinosphere. 



A. The supernova relic and i^e fluxes 

Figure [3] shows the neutrino spectra after propagation in the supernova for the two mass hierarchies and values of the 
third neutrino mixing angle (the magnetic fields or /i^ are set to zero) . The corresponding average neutrino energies 
are shown in Table |lTl Note that the value of 0i3 determines the adiabaticity of the first MSW resonance at high 
density [5l|, H^l , while 9i2 governs the second (adiabatic) one at low density. As expected, when 6*13 is large and the 
hierarchy is normal, maximal conversion occurs for electron neutrinos which acquire the hot spectrum of muon and 
tau neutrinos. On the other hand, for inverted hierarchy (when there is no resonance), the electron neutrino spectra 
are much colder. A similar behaviour is obtained for the electron anti-neutrinos, in the case of the inverted hierarchy. 

The supernova relic Ve and De fluxes, obtained after integrating over the redshift, are shown in Figures S] to [SI The 
results including the RSF effects correspond to /iySo = 0-1 (1) for n = 2 (n = 3). We do not show the supernova 
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NH-L 


NH-S 


IH-L 


IH-S 












MSW 


21.5 17.7 


18.8 17.4 


17.9 21.8 


17.6 18.5 


(^J,B)large 71 = 2 


19.6 20.0 


18.3 20.3 


15.8 18.9 


17.8 14.6 


{^B)large Tl = 3 


19.3 21.4 


18.2 21 


21.5 14.2 


19 16 



TABLE II: Average energies (in MeV) of the Ve and fluxes, after tfie propagation in the SNII environment, with {{^J.B)iarge 
n = 2, 3) EqOor without (MSW) the couphng to the magnetic fields. The different cases correspond to normal (NH) or inverted 
(IH) hierarchy and sin^26li3 = 0.19 (L) or sin^26ii3 = 10"" (S). 





NH-L 


NH-S 


IH-L 


IH-S 






i'e 








Us 




Us 


1 - 50 MeV 


15.7 


10.7 


14.8 


10.7 


15.3 


10.8 


15.3 


10.3 


19.3 - 50 MeV 


1.5 


0.6 


1.2 


0.5 


1.0 


1.0 


1.0 


0.7 


19.3 - 27 MeV 


1.0 


0.4 


0.8 


0.4 


0.7 


0.7 


0.7 


0.5 



TABLE HI: Relic Number Fluxes (cm ^) of electron neutrinos and anti- neutrinos, for different neutrino energy ranges and 
in the MSW case (no interaction with the magnetic field is included). The different cases correspond to normal (NH) or inverted 
(IH) hierarchy and sin^26li3 = 0.19 (L) or sin^26li3 = 10"* (S). 



relic number fluxes obtained for values smaller than 0.001 (0.01) for n — 2 {n — 3) since the RSF effects are small and 
the results are close to the MSW solutions. For values within 0.001 (0.01) and 0.1 (1.) for n = 2 (n = 3) the results 
are in between the RSF and MSW ones. Table Hill presents the relic Ve and De number fluxes integrated over different 
energy intervals. One can see that the relic Ve number flux is larger than the one in practically all cases^. Since 
the relic number fluxes are all peaked in the first 10 MeV, it is clear that this gives the most important contribution 
to the fluxes. Therefore detector energy thresholds around 10 or 20 MeV - as considered or used for example for 
LENA and water Cerenkov detectors respectively - are only sensitive to the tails of the SRN spectra. In the energy 
interval larger than 19.3 MeV, our relic flux ranges between 0.4-0.7 cm"^s~^, depending on the neutrino unknown 
parameters, and is compatible with the Super-Kamiokande limit of 1.2 cm~^s~^ As far as the comparison to 
previous calculations is concerned, our predictions agree well with e.g. [2l|.[l5j. 



B. Events in future observatories 



We consider both the electron anti-neutrino signal associated to scattering on protons and the electron neutrino 
signal due to neutrino capture on nuclei. We take oxygen, carbon and argon as nuclear targets, having in mind the 
MEMPHYS (UNO, Hyper-K), LENA and GLACIER detectors at present under study [3l|. The reactions of interest 
are"^: De +P ^ n + e~^ and fe+^^O-^ e~-|-X or Vg+^^C^ e~-|-"'^^N for the water and scintillator detectors respectively; 
t'e+'^^Ar— > e~-|-'*°K for liquid argon. While the events associated to ^^O cannot be distinguished from the ones on 
protons in water Cerenkov detectors, those corresponding to ^^C can be identified by performing a double coincidence 
with the electrons from ^^N beta decay. It is clear that the both i/e + p ^ n + and the i/e+^"Ar cross sections 
have the advantage of having a very low threshold and of being much larger compared to nuclear targets such as ^^C 
and ^^O (Figure [1]). It is however important to emphasize that while the neutrino-proton cross section is known with 
high accuracy, this is not the case for neutrino-nucleus interactions. One might hope that by the time large-scale 
observatories are in operation, new facilities producing neutrinos in the 100 MeV energy range such as low energy 
beta-beams (60| or those based on conventional sources [6l| might provide us with precise neutrino-nucleus cross 
sections. 

Figures [TKH [9llT0l and fTT]fT2l present the event rates for scattering of De on protons and on oxygen and argon. 



^ Note that our relic Ue fluxes are compatible with the theoretical upper limit of 5.5 cm ■^s ^ extrapolated from the Super-Kamiokande 
result [2I . 

^ Note that we do not consider here events produced by the I'e+^'^C— > e^+^-^B reaction in LENA, since the cross section and the Pe flux 
are smaller than in the fe case. For GLACIER, we do not consider other channels since the charged-current Ue is the dominating one. 
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Neutrino Energy (MeV) Neutrino Energy (MeV) 



FIG. 3: Electron neutrino (left) and anti-neutrino (right) spectra after propagation in a core-collapse supernova, Eqs. ([TJS]) 
with B = 0, for the normal hierarchy with sin^20i3 = 0.19 (NH-L) or sin^26'i3 =10"* (NH-S) (thin solid and dashed lines) and 
the inverted hierarchy with sin^20i3 = 0.19 (IH-L) or sin^ 26li3 =10"* (IH-S) (thick solid and dashed lines). The Fermi-Dirac 
distribution is also shown for comparison (dotted line). 



No Osc. 

— MSWNH-L 

— B: n=2, laige 

— B: n=3, laige 




10 15 20 25 

Neutrino Energy (MeV) 




10 15 20 25 

Neutrino Energy (MeV) 



FIG. 4: Same as Fig. |3]but for the number flux of supernova relic neutrinos and anti-neutrinos Eqs.([T]|8]) with and without 
the coupling to the magnetic fields. Their radial profile is taken as a power law Eq.((5} with n=2,3 and the product fiijB is set 
to 1/0.1 for n = 3/2 (see text). For values much smaller than 0.01 (0.1) for n = 2 (n = 3) the results are close to the MSW 
solution and they are not shown. All curves correspond to the NH-L case. Similar results are obtained for the NH-S. 



respectively. These curves can be easily understood from the the relic number fluxes (Figs. [3] to [S]) and the cross 
sections (Figure [!}, since the neutrino- nucleon and neutrino- nucleus cross sections behave approximately like E^. The 
corresponding expected number of events for De and Ve in 10 years are shown in Tables IIVITW and IVHIVIII for the 
water Cerenkov (MEMPHYS, UNO, Hyper-K), scintillator (LENA^) or hquid argon (GLACIER) technologies. The 
detector sizes are those considered in the corresponding proposals, namely 440 kton fiducial mass (3 shafts), 50 kton 
and 100 kton respectively [3l[. Here we consider that each detector has 100 % efficiency. 

Extensive discussions about possible backgrounds can be found in the literature. A detailed investigation is pre- 
formed in 21] for water Cerenkov, in [TtI for scintillators and in [26] for liquid argon. This issue is clearly crucial 
for the observation of relic supernova neutrinos. Three are the main sources of the backgrounds: solar neutrinos, 



* Note that the Ue+p event rates presented here for LENA are obtained with = 1 while in ilTII with = 2.5. 

^ LENA is composed by 20 % PXE and 80 % Dodecane. The number of protons in the fiducial volume is 2.9 X 10''^. 
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FIG. 5: Same as Fig. |4]but for electron neutrinos in the IH-L (left) and IH-S (right) cases. 
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FIG. 6: Same as Fig. [5] but for electron anti-neutrinos. 



reactor anti-neutrinos, atmospheric neutrinos. Solar neutrinos have low energies, the maximal being of almost 20 
MeV coming from hep. Since water Cerenkov detectors do not have any charge identification, electron neutrinos can 
be mis-identified as anti-neutrinos. However, this background can be almost fully suppressed thanks to the good 
angular resolution of such detectors. On the other hand, the signal from electron neutrinos and anti-neutrinos can be 
distinguished in scintillator detectors. Concerning reactor electron anti-neutrinos, the flux is peaked at 2-3 MeV and 
has a tail reaching about 13 MeV. This contribution might be important or not depending on the chosen location. If 
present, it can clearly hinder SRN coming from z >l-2. In the case of the capture of the relic electron neutrinos by 
oxygen or carbon, the reaction threshold is so high that neither reactor nor solar neutrinos can produce a background. 
Atmospheric neutrinos produce events at higher energy. Note that in water Cerenkov detectors there is a supplemen- 
tary source of background due to invisible muons. These are produced by atmospheric neutrinos that interact with 
the nucleons in the detector producing muons with energies below the Cerenkov radiation threshold. Figure [T51 shows 
as an example the backgrounds in a scintillator detector such as LENA in comparison with our 9^ rates. 

Various authors have identified an energy window for the observation of relic supernova neutrinos, where the 
backgrounds are lower than the expected signal. At the present stage it might be premature to define an energy 
window for two reasons. First some backgrounds like for example the atmospheric one might scale differently in large- 
scale observatories from what is currently assumed. Secondly there are still many theoretical uncertainties coming 
either from supernova modelling, or from the knowledge of the supernova rate and from still unknown neutrino 
properties, in particular the knowledge of the hierarchy or of the ^13 value. For the above mentioned reasons we 
present the expected rates in two different energy windows. The 6.8 and 5-50 MeV energy ranges^are taken to show 



The threshold of 6.8 MeV takes into account the 1.8 MeV threshold energy for the reaction Ue+p — >n + e+ plus the 5 MeV necessary 
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Neutrino Energy (MeV) Neutrino Energy (MeV) 



FIG. 7: Event rates (100 kton x year)~^ for + p — > n + e"*" due to SRN for a normal hierachy and a large (left) or small 
(right) value of 613. The various curves present the results when the interaction with matter (MSW) and the coupling of the 
neutrino magnetic moment to the stellar magnetic fields are included (B large) with a radial dependence of r~^ (n=2) or r~^ 
(n=3) Eq.©. The expected events when no oscillation is taken into account are shown for comparison. 




Neutrino Energy (MeV) Neutrino Energy (MeV) 



FIG. 8: Same as Figure [7] but for the inverted hierarchy. 



how the event rate might increase if one manages to have a threshold as low as 6.8 and 5 MeV for the outgoing lepton 
and as high as 50 MeV, by suppressing the atmospheric background. Such an information, that might be helpful in the 
future, should be taken as indicative only. The 19.3-27, 9.3-25, 17.5-41.5 MeV for water Cerenkov, scintillators, liquid 
argon detectors are taken following Refs. [2l[ , [13] , [11] respectively. We make this choice to facilitate the comparison 
with available predictions^. 

As can be seen from Tables llVlfVl in absence of the coupling to magnetic fields the largest number of events is 
expected for the IH-L case because of the MSW-/i resonance. The events are larger by 50-60 %. In the normal 
hierarchy (^13 large or small) case, electron anti-neutrinos undergoes regular matter effects and the RSF-Zi resonance 



to observe the outgoing electron in water Cerenkov. The threshold of 5 MeV for scintillators is fixed in the same way. 
^ Note that our relic fluxes cross the backgrounds at higher energy (see for example I13t . This would give energy ranges larger than the 
presented ones. 



10 



0,03 




0,03 



20 30 40 

Neutrino Energy (MeV) 



50 



0,01 



No Osc. 
MSW NH-S 
B: n=2, large 
B: n=3, large 




20 30 40 

Neutrino Energy (MeV) 



50 



FIG. 9: Same as Figure[7]but for the event rates (100 kton x year) ^ associated to Ve + ^^O^ e + X. 
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FIG. 10: Same as Figure [9] but for the inverted hierarchy. 



which renders the electron anti-neutrinos hotter compared to the case with matters effects only. Therefore in presence 
of strong magnetic fields and/or large values of the event rates can be increased by 30 to 50 %. For the inverted 
hierarchy, electron anti-neutrinos encounter both the MSW-/i and the RSF-/i resonances. In this case the events are 
significantly reduced. The two magnetic field profiles considered give similar results, except for the inverted hierarchy 
and large mixing angle. For this case and if the n — 2 profile is taken, the anti-neutrinos undergo the RSF-/ resonance 
as well. As a consequence the number of events is doubled compared to the n = 3 profile and is close to the MSW 
result. 

In Tables IVHIVIII we show the predicted events in the case of electron neutrino-nucleus scattering, choosing large 
intervals only^, except for the argon case. In fact, this is the first time that interactions of SRN on oxygen and carbon 
are considered. Since a detailed study of the atmospheric backgrounds for these cases does not exist, we hope that our 
work will trigger investigations of this channels for MEMPHYS (UNO, Hyper-K) and LENA. A few events are found 
for z^e+^^C, up to 20 for and 100-200 for 1/6+"*° Ar. Note that the large number of events found for 1^6+"*° Ar 



* The 20.5 MeV threshold includes the reaction threshold and 5 (3) MeV electron energy for MEMPHYS (LENA). 
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FIG. 11: Same as [7] but for the event rates (10 kton x year) ^ associated to Ue+'^^' At~* e +^''K 
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FIG. 12: Same as Figure [TT] but for the inverted hierarchy. 



is due both to the large cross section (Figure [T]) and to the large i^e flux (Table HII)) . In absence of magnetic moment 
effects, the maximal conversion for electron neutrinos occurs for a normal hierarchy and a large 6*13 value, due to the 
MSW-/i resonance. In fact, since gets a hotter spectrum because of the mixing with v^, neutrinos the number of 
events is increased by 30-60 % compared to the case when the resonance is not (fully) adiabatic or of normal matter 
effects. On the other hand, RSF conversions from j/g to i?^, Vr generally decrease by 20-40 %, compared to the MSW 
only, in most cases. The only exception is for the inverted hierarchy and the n = 3 magnetic field profile for which 
the events increase by 15 up to 50 %. It is clear that the search for relic might either improve the present LSD 
limit [5*] (in case of non-observation) or bring crucial information to our understanding of core-collapse supernova and 
neutrino physics. 
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FIG. 13: Expected event rates in the LENA detector for normal (NH) and inverted (IH) liierarcliy as well as a small (S) and a 
large (L) third neutrino mixing angle. Both the reactor electron anti-neutrinos (dotted line) at Pyhasalmi and the atmospheric 
(dashed-dotted line) backgrounds are shown. The backgrounds are from Ref.[l3,|62. 



MEMPHYS LENA 




6.8 - 50 


19.3 - 27 


5.0 - 50 


9.3 - 25 


MSW NH-L 

[^J.B)large NH-L 


652 
873 (968) 


142 
214 (240) 


67 
89 (97) 


44 
55 (61) 


MSW NH-S 

itlB)larse NH-S 


631 
868 (921) 


135 
209 (225) 


65 
88 (94) 


43 
55 (58) 



TABLE IV: Expected number of events after 10 years for electron anti-neutrino scattering on protons in MEMPHYS (water 
Cerenkov) or LENA (scintillator) detectors. The size considered is 440 kton fiducial mass, 50 kton and 100 kton respectively, 
following the corresponding proposals [3l[. An efficiency of 100 % is assumed. The calculations correspond to the normal 
hierarchy (NH) and large (L) or small (S) value of the third neutrino mixing angle. The calculations include the coupling 
of neutrinos to matter effects without (MSW) or with the coupling to the stellar magnetic fields ((/i-B)!arge) . The stellar 
magnetic profile is taken to be a power- law Eq.([5ll with n = 2 (n = 3 results are given in parenthesis). Two energy intervals 
are considered (see text). 



MEMPHYS LENA 




6.8 - 50 


19.3 - 27 


5.0-50 


9.3-25 


MSW IH-L 

{lJ-B)targe IH-L 


997 
776 (464) 


248 
181 (95) 


101 

80 (48) 


64 
51 (30) 


MSW IH-S 

{)J,B)large IH-S 


751 
479 (559) 


173 
98 (120) 


77 
51 (58) 


50 
32 (38) 



TABLE V: Same as Table |IV] but for the inverted hierarchy. 



IV. CONCLUSIONS 



In this work we have analysed both the relic supernova electron anti-neutrino and neutrino fluxes and the associated 
signals due to scattering on protons and nuclei respectively, in large-scale observatories such as MEMPHYS (UNO, 
Hyper-K), LENA or GLACIER. Our calculations of the neutrino propagation in the star include both the matter 
interaction and the resonant spin-flavour conversion. One of our goals has been to try to investigate under which 
conditions (profile, fiij and/or Bq values) RSF effects could be seen. In all cases studied, the results show that RSF 
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TABLE VI: Expected number of events after 10 years for electron neutrino scattering on nuclei in MEMPHYS (water Cerenkov), 
LENA (scintillator) or GLACIER (liquid argon) detectors. The size considered is 440 kton, 50 kton and 100 kton respectively, 
following the corresponding proposals 31]. An efficiency of 100 % is assumed. The calculations correspond to the normal 
hierarchy (NH) and large (L) or small (S) value of the third neutrino mixing angle. The calculations include the coupling 
of neutrinos to matter effects without (MSW) or with the coupling to the stellar magnetic fields {{fJ,B)targe) ■ The stellar 
magnetic profile is taken to be a power-law Eq.([5)l with n = 2 (n = 3 results are given in parenthesis). Two energy intervals 
are considered in the case of argon (see text). 
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10 (14) 




1.5 (2.1) 




181 (221) 


112(142) 



TABLE VII: Same as Table |Vll but for the inverted hierarchy. 



phenomenon might increase/decrease significantly the number of events, for strong magnetic fields and/or large values 
of the neutrino magnetic moment. It is interesting to note that these effects are maximal for electron anti-neutrinos, 
when no resonant conversion is present due to matter efffects (normal hierarchy and a large or small value of ^13) 
due to the RSF-ft. resonance. In the case of electron neutrinos, the RSF tend to decrease the events, in most cases. 
Our results are also sensitive to the magnetic field profiles. To obtain thes results with the spin-flavour conversion, 
we have made the assumption of a static density supernova profile and also neglected possible effects of the shock 
wave on the stellar magnetic fields. Our present knowledge in this respect is in fact still too limited to include such 
aspects. From this point of view, our analysis should be taken as a first step in this direction and should be updated 
once progress is performed on the knowledge of the supernova magnetic fields and on the shock wave impact. 

Our main goal has been to perform a numerical study of the relic fluxes and their signal in the next generation 
observatories. We have obtained that z/g fluxes are larger than the De up to a factor 2.5 depending on the oscillation 
scenarios (normal or inverted hierarchy and ^13 value). We have shown that, after 10 years operation, between 
100 and 250 i>e + P events are expected in the 19.3-27 MeV neutrino energy range for MEMPHYS, while 45-65 are 
expected in LENA in the 9.3-25 MeV energy range. We flnd that neutrino-carbon scattering will produce 2-3 events, 
while neutrino-oxygen will give 10-20 events and the neutrino-argon interaction will induce 100 to 200 events. (These 
number are obtained by taking the normalization factor for the star formation rate = 1, which implies a local 
formation rate of 0.7x lO^^/iyoM© year~^Mpc~'^, consistent with mildly dust-corrected UV data at low redshift). The 
possible measurement of both z^e and De could help us determining the neutrino fluxes at the neutrinosphere and/or 
still unknown neutrino properties, such as the hierarchy or the knowledge of the third neutrino mixing angle. This 
will depend on the progress made in supernova modelling or neutrino physics by the time such observatories might 
come into operation. In conclusion, future large-scale observatories based on water, scintillator or liquid argon have 
a great discovery potential for relic supernova neutrinos. Their observation would clearly provide us with essential 
information on core-collapse supernova explosions, the supernova explosion rate and important neutrino properties 
that still remain unknown. 
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